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Abstract. We report on a detailed analysis of a SO(IO) SUSY GUT model of Dermfsek and Raby 
(DR) with a D3 family symmetry. The model is completely specified in terms of only 24 parameters 
and is able to successfully describe both quark and lepton masses and mixings, except for \Vub\ 
that turns out to be too low. However, a global fit shows that flavor changing (FC) processes like 
Bs — *■ ijf^ n~ , Bs-mixing, — > t^v, B — > Xs7 and B Xs£'^£~ pose a serious problem to the 
DR model. The simultaneous description of these FC processes forces squarks to have masses well 
above 1 TeV, not appealing on grounds of naturalness and probably beyond the reach of the LHC. 
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1 Introduction 

Many extensions of the Standard Model (SM), like 
the general Minimal Supersymmetric Standard Model 
(MSSM), typically introduce a large set of additional 
parameters to the SM ones and therefore largely lose 
their predictivity. On the other hand, in Supersymmet- 
ric Grand Unified Theories (SUSY GUTs) with addi- 
tional family symmetries, the number of parameters 
can even be smaller than in the SM. In a top down ap- 
proach, these models then allow to predict observables 
at the low scale in terms of a manageable number of 
GUT scale parameters. 

One of such highly predictive models is the SO(IO) 
SUSY GUT introduced by Dermfsek and Raby in 
In [2] lepton flavor violating processes and electric di- 
pole moments were studied extensively within this mo- 
del. Here we report on a detailed analysis [3j of the 
same model in the light of the best measured FC pro- 
cesses in the quark sector. 



2 The Model 

The DR model is a supersymmetric SO(IO) Grand 
Unified Theory that is supplemented by a D3 x[U{l)x 
Zi X Z3] family symmetry. 

The three generations of quarks and leptons are 
each unified in a 16. Furthermore, the model has one 
additional 10 that contains the two Higgs doublets of 
the MSSM. The family symmetry ensures that only a 
universal third generation Yukawa coupling is allowed. 
Yukawa couplings for the first and second generation 
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Table 1. The 24 parameters in the DR model. (The 
Yukawa textures p, a, e and ^ are complex.) 



Sector 


# 


Parameters 


Yukawa textures 


11 


e, e', A, p, a, e, i 


RH neutrinos 


3 


Mr,, Mr^, Mr^ 


gauge couplings 


3 


aa, Ma, es 


SUSY (GUT scale) 


5 


M-i/2, mi6, Ao, niH^, mn^ 


SUSY (EW scale) 


2 


tan /3, /i 



are then generated by a Froggatt-Nielsen mechanism 
4 . It turns out that the resulting Yukawa matrices can 
be parameterized in terms of only 11 parameters. All 
parameters of the model are summarized in Table [TJ 
Among them, the most important role in the numerical 
analysis of the FC processes is played by the universal 
sfermion mass mig, the universal trilinear coupling Aq 
and the Higgsino mass parameter /i. Moreover, tan/3 
is forced to be around 50, because of third generation 
Yukawa unification. 

The total number of model parameters is 24 and 
once they are fixed, the complete MSSM Lagrangian 
at the electro- weak scale is specified. 

3 Basic Procedure of the Analysis 

Starting with the model parameters at the GUT scale, 
the Yukawa matrices, the right-handed (RH) neutrino 
mass matrix, the gauge couplings and the soft SUSY 
breaking parameters are run down using renormaliza- 
tion group equations. The RH neutrinos are integrated 
out at their respective scale and the remaining param- 
eters are further run down to the electro-weak scale. 
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Fig. 1. Schematic cliart of the strategy followed in the 
fitting procedure of [3]. 



where also the parameters ^ and tan/3 are specified. 
The complete set of MSSM parameters is then given in 
terms of the original model parameters from Table [T] 
This is especially true for the flavor off-diagonal en- 
tries of the squark mass matrices that are generated 
radiatively in the running procedure because of the 
Yukawa couplings. As the model features no additional 
sources of flavor violation and in particular no other 
CP phases than the ones appearing in the CKM and 
PMNS matrices, it can be classified as being minimal 
flavor violating [5]. 

Having at hand all MSSM parameters, one calcu- 
lates the SUSY spectrum, loop corrections to fermion 
masses and mixings and finally also FCNC observ- 
ables like AMs, BR(Bs ^ /U+M"), BR(B Xs-y) and 
BK{B Xst^i^) as well as the branching ratio for 
the FC decay B+ r+i/. 

Using these FC observables as well as flavor con- 
serving quantities like gauge couplings and fermion 
masses and mixing angles, a function is defined [3], 
that is then minimized by varying the original model 
parameters. (See Fig. [T] for a schematic chart of this 
whole procedure.) 

Using only fiavor conserving observables, it was 
shown in [I] that the DR model is indeed able to give 
successfull fits to the experimental data. In particular, 
the model is able to excellently reproduce the CKM 
matrix elements except for Vub, whose absolute value 
typically comes out as 



IK 



DRi 



3.2 X 10 



-3 



(1) 



This number is smaller than both the value from the 
exclusive and the inclusive determination [6] 



I Kb 



(3.50 ±0.40) X 10" 



K;,|i„Ti= (4-49 ±0.33) x 10" 



(2) 



The main novelty of the analysis in [3] was then to 
also include the abovementioned FC observables into 
the fit. 



4 Interplay of Flavor Changing Processes 

Although the DR model is minimal flavor violating, 
one expects interesting effects in various FC processes 
due to the large value of tan/jQ In this section we 
discuss the general pattern of these effects. 



4.1 Bs — > fJ-'^fJ' and Bg mixing 

Combining data from CDF and D0 results in the fol- 
lowing upper bound on the branching ratio of the rare 
decay Bg M^A*^ 95% C.L. 

BR{Bs M+A'")"'"' < 5.8 x 10-^ , (3) 
that is still much larger than the SM prediction [5113] 
BR{B, M+/i")^^ = (3.37 ± 0.31) x 10"^ . (4) 



The helicity suppression of the SM result can be lifted 
in the MSSM with large tan/3 by neutral Higgs pen- 
guins [9] that lead to contributions to the branching 
ratio that are strongly enhanced by tan /3 



tan6/3 



(5) 



The same mechanism also leads to neutral Higgs dou- 
ble penguin contributions to the mass difference in the 
Bs — Bs system [7] 



(AMs) 



DP 



tan4/3 



(6) 



On the experimental side, this quantity is known very 
precisely [TU] 



(AMsY^'P = (17.77 ± 0.10 ± 0.07)ps~i 



(7) 



On the other hand, the theory prediction in the SM 
suffers from large hadronic uncertainties l6] 



(AMsf^ = (18.6±2.3)ps-i , 



(8) 



leaving still some room for new physics contributions. 
But as tan /3 is forced to be around 50 by third gener- 
ation Yukawa unification, both observables constrain 
the pseudoscalar Higgs mass Ma- In fact, in the DR 
model we find a lower bound on > 450 GeV, 
that then approximately also holds for the other heavy 
Higgs particles. 



4.2 B+ T+u 

Using the most recent experimental results one obtains 
the following average for the branching ratio of the tree 
level decay 5+ t^v (see and references therein) 



BR(B^ 



cxp 



(1.41 ±0.43) X 10" 



(9) 



^ In the numerical analysis, we resum large tan/3 correc- 
tions following l^. 
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The SM branching ratio is proportional to |T46p- Us- 
ing the exclusive and the inclusive value for | Vub \ from 
eq. ([2]) yields the following SM predictions [3j 



BR(B+ 
BR(B+ 



(0.80 ±0.20) X 10-'' 
(1.31 ±0.23) X lO"'' 



(10) 



In the MSSM there is an additional contribution to this 
decay coming from the exchange of a charged Higgs 
boson. It interferes destructively with the SM contri- 
bution and one finds [TT 1 I12 1 I3] 



Rbtu 



BR(B^ 



BR(B^ 



SM 



= 1- 



M2 



tan^/? 



Ml, 1 



£0 tan P 



* uh 



T/SM 



(11) 



As confirmed in ^ , the value for \Vub \ in the DR model 
is always lower than in the SM, which leads to a further 
suppression of the branching ratio with respect to the 
SM value. In the DR model, we typically find 



BR(B^ 



< 0.6 X 10-4 



(12) 



which is however not yet excluded, given the large ex- 
perimental error in eq. ([9|). 



4.3 B and B Xsl+l' 

The experimental value for the branching ratio of the 
inclusive decay B — > reads [13j 



BR(B ^ Xsi) 



cxp 



(3.55 ±0.27) X 10" 



(13) 



which is sUghtly above the NNLO SM prediction [Mj 



BR(B ^ Xsi) 



SM 



(3.15 ±0.23) X 10" 



(14) 



As discussed in Sec. 14.11 Higgs masses are forced to be 
quite large, implying that new physics contributions 
to C7, the Wilson coefficient governing B Xgj, are 
dominated by chargino - stop loops. For large values of 
tan P these chargino contributions obey the following 
approximate relation [T51IT2] 



cx fiAt tan /3 x sign(Cf'^') 



(15) 



For /i > and At < the sign of the chargino contribu- 
tion is opposite to the SM one. Without invoking fur- 
ther constraints, the model favors very large chargino 

contributions C^"" « -2C^^ that lead to C7 « -C^^ 
which accommodates the data on _B — s- X^^. 

A further important process to be considered is 
then B Xs(~^£~- Both the forward backward asym- 
metry and the branching ratio of this decay are sensi- 
tive to the sign of C\. In case the sign of C7 is oppo- 
site to its SM value, the forward backward asymme- 
try has no zero, which is however not yet excluded 
experimentally. On the other hand, the experimen- 
tal data on the branching ratio in the low s region 
1 GeV^ < s < 6 GeV^ Ql 



BR(B ^ Xsi+e-y^'P = (1.60 ± 0.51) 



10" 



(16) 
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Fig. 2. Total vs. the lightest stop mass for all obtained 
fits. Red circular points correspond to fits with positive /i, 
blue squares to negative /j.. 



is in very good agreement with the SM prediction [17j 

BR(B ^ Xsi+r f^ = (1.59 ± 0.11) X 10"*^ . (17) 

It has been shown p^SlfTS] that the experimental result 
(fTB)) excludes the "wrong sign" solution for C7 if the 
Wilson coefficients Cg and Cio are SM-like. This is es- 
pecially the case in a minimal flavor violating MSSM 
[l9] and also in the DR model. Thus chargino contri- 
butions to C7 have to be suppressed, which can only 
be done by raising the stop masses. 



5 Results of the Numerical Analysis 

The main features of the interplay of the FC processes 
described in Sec. [4] are then also reflected in the per- 
formed numerical fits. The adopted strategy in [3] was 
to roughly set the scale for the sfermion masses by fix- 
ing mi6, while all the other model parameters where 
left free in the fits. 



5.1 Fits with fj. > 

For positive values of /i, the fit strongly prefers values 
for Aq that obey the following approximate relation at 
the GUT scale 

Ao « -2mi6 , (18) 

which helps to obtain third generation Yukawa uni- 
fication [301 and leads to an inverted mass hierarchy 
for squarks. These large values for Aq also result in 
large negative values for At at the electro- weak scale, 
that in turn lead to the large chargino contributions 
to i? ^ Xsj discussed in Sec. 14.31 The only possibihty 
to tame these corrections is then to decouple stops, 
which can be done by choosing a very large mig . 
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Table 2. Fit results for some selected observables for two 
obtained fits. More detailed tables can be found in [3J. 

mi6 = 10 TeV, = 1.2 TeV 



Observable 


Fit value 


Pull (a) 


mi, [TeV] 


1.9 




BR(B, M+M") X 10* 


2.1 




BR(B+ ^ T+u) X 10* 


0.517 


2.1 


BR(B ^ Xs-y) X 10* 


2.86 


1.3 


mi6 = 4 TeV, /i = -2.1 TeV 


Observable 


Fit value 


Pull (a) 


mi, [TeV] 


2.6 




BR(B, /i+^i") X 10* 


0.33 




BR(B+ r+zy) X 10* 


0.59 


1.9 


BR(B -> X,7) X 10* 


3.34 


0.4 



In Fig. [2] the correlation between the lightest stop 
mass rrif^ and the total is shown for all obtained fits. 
Reasonable fits require stop masses TOj^ > 1.9 TeV, 
corresponding to values for mig in the range (8 — 
10) TeV. For such heavy stops, the chargino correc- 
tions to -B — > Xs"f are under control, but the large 
stop masses clash with the motivation for SUSY as a 
solution of the hierarchy problem. 

The fit results for some observables in this scenario 
are collected in Table [21 



5.2 Fits with < 

In [3] also scenarios with negative were considered. 
In such cases Aq and mie apparently do not have to 
fullfill relation p^ PI The fit then chooses values for 
Aq that typically lead to very small At- This makes 
chargino contributions to C7 automatically small (see 
eq. (|15p ). thus solving the problem with B — > Xg'y. On 
the other hand, the squark spectrum does not show an 
inverted hierarchy in this case and the lightest squark 
(which is usually still a stop) has again a very large 
mass m^^ > 1.8 TeV. 

Some fit results for a negative /i case can again be 
found in Tabled 



6 Conclusions 

The SO(IO) SUSY GUT model of Dermisek and Raby 
[1] is able to successfully fit the known quark and lep- 
ton masses as well as the CKM and PMNS mixing 
matrices. The only exception is the absolute value of 

^ Successfull fits away from relation (I18|l were obtained 
also for positive values of /i. However, these fits perform 
worse than the corresponding negative fj, fits, because large 
loop corrections to the bottom quark mass add up con- 
structively and lead to a prediction of mt that is roughly 
4(T too large. 



the CKM matrix element Vub that is even smaller than 
the central exclusive value. 

Given such a small value of \Vub\, we then find a 
very low upper bound (jl2p on the branching ratio of 
^ T+i/ in the DR model. Consequently, this decay 
will turn out to be quite problematic for the model, if 
the central experimental value for the branching ratio 
stays above 1.0 x 10"''. 

Furthermore, we find that the model is not able to 
simultaneously fit the branching ratios of the decays 
Bs fJ.~^t^~, B and B XJ+£-, unless the 

squark spectrum is made very heavy (m^^ > 1.8 TeV). 
Such large squark masses may be problematic from the 
point of view of naturalness and squarks may be even 
beyond the reach of the LHC. 

As the example of the DR model shows, it is es- 
sential to check simultaneously many flavor changing 
processes to test the validity of models for fermion 
masses and mixings. 
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